A model of HIV budding and self-assembly, role of cell membrane 



OO 

o 
o 

(N 



u 

6 

x> : 



> 
o 

00 
(N 

cn 

o 

OO 
O 



X 



Rui Zhang and Toan T. Nguyen 
School of Physics, Georgia Institute of Technology, 837 State Street, Atlanta, Georgia 30332-0430 

(Dated: March 20, 2008) 

Budding from the plasma membrane of the host cell is an indispensable step in the life cycle 
of the Human Immunodeficiency Virus (HIV), which belongs to a large family of enveloped RNA 
viruses, retroviruses. Unlike regular enveloped viruses, retrovirus budding happens concurrently 
with the self-assembly of retrovirus protein subunits (Gags) into spherical virus capsids on the cell 
membrane. Led by this unique budding and assembly mechanism, we study the free energy profile of 
retrovirus budding, taking into account of the Gag-Gag attraction energy and the membrane elastic 
energy. We find that if the Gag-Gag attraction is strong, budding always proceeds to completion. 
During early stage of budding, the zenith angle of partial budded capsids, a, increases with time 
as a cx t 1 / 3 . However, when Gag-Gag attraction is weak, a metastable state of partial budding 
appears. The zenith angle of these partially spherical capsids is given by oiq ~ (t 2 /k<t) 1 / 4 in 
a linear approximation, where k and o are the bending modulus and the surface tension of the 
membrane, and r is a line tension of the capsid proportional to the strength of Gag-Gag attraction. 
Numerically, we find qo < 0.3-7T without any approximations. Using experimental parameters, we 
show that HIV budding and assembly always proceed to completion in normal biological conditions. 
On the other hand, by changing Gag-Gag interaction strength or membrane rigidity, it is relatively 
easy to tune it back and forth between complete budding and partial budding. Our model agrees 
reasonably well with experiments observing partial budding of retroviruses including HIV. 



I. INTRODUCTION 

The Human Immunodeficiency Virus (HIV) is famous 
for its ability to induce Acquired Immunodeficiency Syn- 
drome (AIDS). It belongs to a large family of enveloped 
RNA viruses, retroviruses. Retroviruses are character- 
ized by the unique infection strategy of reverse transcrip- 
tion, in which the genetic information flows from RNA 
back to DNA (therefore the name "retro" ) [l[ . Budding 
is an indispensable step in the retrovir al life cycle @,||. 
After the major retroviral structural protein, Gags, are 
synthesized inside the host cell, they are transported to 
the cell membrane and self-assemble into spherical pro- 
tein shells called "capsids" , with viral RNA genome and 
other auxiliary viral proteins packaged inside. At the 
same time, these capsids, enveloped by the cellular mem- 
brane, have to bud out of the membrane to target other 
host cells. In other words, budding and assembly of 
retroviruses happen concurrently on the cell membrane. 
Despite a large body of experiments done within the last 
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FIG. 1: Electron microscopic image of partial budding of 
HIV-1 viruses (reprinted from [2J). 



decade, the biological pathway and mechanism of retro- 
viral budding have still not been fully understood 0, 0| • 
One important unexplained observation is that viral bud- 
ding can be inhibited partially or completely by chang- 
ing the cell environment or mutating the late domains of 
the Gag proteins. In these situations, capsids are only 
partially formed and stuck on the membrane (Fig. [T]). 
Motivated directly by this partial budding phenomenon, 
in this paper, we propose a physical model to study HIV 
(and retroviruses in general) budding and assembly on 
the elastic membrane. Physically, this situation is in- 
teresting because it provides a unique two dimensional 
self-assembly mechanism in which the membrane elastic 
energy plays an important role, since assembly is always 
accompanied with budding. Biologically, understanding 
the physical mechanism of HIV budding and assembly 
is certainly important toward understanding the HIV 
life cycle. It is also important in the light of recent ef- 
fort from the virology community to develop assembly- 
oriented anti-viral therapy. 

Budding of regular enveloped viruses was studied the- 
oretically by several authors (TDGB) in Ref. 0, i, i]. 
However, the viral budding pathway and the physical 
model studied by TDGB is qualitatively different from 
retroviral budding we study in this paper. For regular 
enveloped viruses, viral capsids are fully assembled in- 
side the cell 0,11, Q- After that, they are transported to 
the cell membrane, bind to the viral spike proteins (em- 
bedded in the cell membrane) and then bud out through 
the cell membrane (see Fig. Od). Therefore, capsids for- 
mation and budding out of the membrane are two sepa- 
rate processes. And the main driving force of budding is 
the capsid-membrane attraction (mediated by embedded 
spike proteins). 

Budding of retroviruses follows a completely different 
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FIG. 2: (Color online) Schematic illustrations of two different 
types of virus budding, (a): budding of retroviruses: cap- 
sid proteins (Gags) are first attracted to the membrane, then 
self-assemble and bud on the membrane at the same time. 
This is the system studied in this paper, (b): budding of 
regular enveloped viruses: capsid proteins first self-assemble 
into complete capsids inside the cell, then bud on the mem- 
brane Q. (a) also shows the cylindrical coordinate system 
(r, h, </)) used in the model (the polar angle <j> is not shown). 



pathway. Various TEM and X-ray tomography experi- 
ments suggest that retroviral capsids are assembled from 
Gag proteins on the cell membrane and bud out of the 
cell concurrently @, Q- Based on these experiments, we 
study a different model for HIV (and retrovirus in gen- 
eral) budding and assembly shown in Fig. In this 
new model, we assume retroviral capsids are assembled 
from membrane-bound Gags only, neglecting the possi- 
bility that Gags from the interior of the cell may par- 
ticipate. In other words, the Gag-membrane attraction 
is strong such that Gags always bind to the membrane. 
This assumption is supported by various experimental 
observation where budding is completely inhibited (no 
capsids are formed ) bu t Gags are found in abundance at 
the cell membrane [l0| . In contrast to the TDGB model, 
the primary driving force of our retroviral budding is the 
short range attraction between these membrane-bound 
Gag proteins. This correlates well with experimental fact 
that point mutations changing Gag-Gag interactions af- 
fect the degree of viral budding. On the other hands, 
spike proteins or virus RNA seem not important for 
retroviral budding. In vitro, Gag proteins are directly 
attracted to the membrane and they alone are usually 
sufficient for the assembly and release of virus-like parti- 
cles [1, [ll|, [lH . We therefore neglect the contribution of 
all other proteins or RNA components of retroviruses in 
our model. 

In this paper, for a given set of parameters (the mem- 



brane Gag concentration, the Gag-Gag interaction, and 
the cell membrane bending and stretching rigidity), we 
study the free energy profile of budded viral capsids. 
Two energies are considered explicitly: first, the elas- 
tic energy of the membrane including the bending and 
stretching energy; second, the Gag-Gag attraction en- 
ergy when a Gag makes contact to the other Gag (see 
Fig. [2^). Since the elastic energy scale is much larger 
than ksT, for example, the bending rigidity of normal 
membranes is about 20 keT, thermal fluctuations are 
higher order corrections and neglected in the theoretical 
treatment. Focusing on the budding process, we also as- 
sume that the Gag-Gag interaction is strong enough such 
that the entropic cost of bringing free Gags to the capsid 
can be ignored. For simplicity, we assume the shape of 
the capsid together with the membrane attached to it is 
(partially) spherical with radius R (Fig. \2k)- The size 
of a capsid is then characterized by the zenith angle a 
at its edge, the smallest being the angle of a single Gag 
protein, ao (Fig. [5^). Since ac is very small (cvg = 0.03 
for a typical HIV capsid containing 5000 Gags), we take 
ag — > in the theoretical consideration and treat a as 
a continuous variable. As budding proceeds to comple- 
tion, a increases from ac to ir. When a — n, the capsid 
actually leaves the membrane through membrane fission. 
In this paper, we do not consider this fission process and 
thus, in our terminology, complete budding always means 
a — ► it. To simplify the calculation, we employ an scaling 
description where we neglect the variation in the degree 
of viral budding and assume all capsids have the same 
average zenith angle a. 
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FIG. 3: The schematic illustration of the total free energy 
density as a function of the capsid size a. The left and right 
profiles correspond to strong or weak Gag-Gag attraction re- 
spectively. Here r is the line tension of the rim of a partially 
budded capsid. r is proportional to the strength of Gag-Gag 
attraction. r c is the threshold line tension at which the local 
minimum at Qo appears. 

Our main result is shown in Fig. [3] The key parameter 
is the strength of Gag-Gag attraction which can be ad- 
justed experimentally by mutating of Gags, complexing 
of Gags with other molecules or by changing pH or salin- 
ity of the cell cytoplasm near the membrane [3, fl3j . In a 
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partially budded capsid, the line tension r of the rim of 
a capsid is directly proportional to this interaction Gag- 
gag interaction. When the Gag-Gag attraction is strong 
(or when r is greater than a threshold value r c ), as in 
the normal biological conditions of HIV, budding always 
proceeds to completion, i.e., a — > tt (the left panel of 
Fig. [3]) . At early stage of budding, the size of a partially 
budded capsid increases very slowly with time: 



«(*) « Wr d iff) 1/3 



(1) 




where the time scale Tdi / / depends on the lateral mobility 
of Gag, the radius of the capsid and the initial concen- 
tration Gag (see Eq. I|5ip). On the other hand, when the 
Gag-Gag attraction is weak (r < r c ), for example, after 
mutation of the late domains of the Gag protein, partial 
budding appears as a metastable state at the capsid size 
«o (the right panel of Fig. [3]). In this case, the free en- 
ergy barrier can be much larger than fc^T and budding is 
kinetically trapped at oeo. Using a linear approximation, 
we find 



(2) 



where a and k are the surface tension and bending rigid- 
ity of the membrane. 

The energetics of HIV budding and assembly is studied 
both analytically and numerically in this paper. Analyt- 
ically, the complete scaling behaviors of the free energy 
density profile in asymptotic limits of "soft" and "stiff" 
membranes are calculated (the meaning of "soft" and 
"stiff" membrane will be clear in later sections). In all 
cases, they agree with the numerical result well. How- 
ever, the numerical result gives a complete solution to 
the problem including nonlinear regimes where the an- 
alytical result is normally not available. The inequility 
ao < 0.37T is found to always hold from the numerical 
calculation without any approximations. 

It is worth to point out that budding in our model can 
be considered as a consequence of the inhomogeneity of 
the membrane if one considers Gags as a part of the mem- 
brane. In this sense, our work is related to Julicher and 
Lipowsky's works on domain-induced budding of vesi- 
cles [T3.ll5|. However, in their papers, the inhomogeneity 
was introduced through two kinds of lipids which do not 
carry a given curvature like our Gags. Domain-induced 
budding is a consequence of demixing of these different 
molecules. As a result, their budding happens in a much 
larger length scale (comparable to the size of the vesi- 
cle) where only two phases coexist, one budded out from 
the other. While in our case, wc consider budding at a 
much smaller length scale (a typical HIV-1 virus particle 
is about 130 nm in diameter, which is a hundred times 
smaller than the size of a host cell) and actually have 
a multi-phase coexistence since there are more than one 
capsid on the membrane. 

This paper is organized as follows. In Sec.|TT], we intro- 
duce the physical model of HIV budding and assembly. 



We then discuss the analytical solution to the elastic en- 
ergy of the membrane in Sec. IIIII and to the total free 
energy density in Sec. IIVI The numerical result is then 
provided and compared to the analytical results in Sec.lVl 
After we get the complete theoretical result, we discuss 
budding kinetics and make connections to experiments in 
Sec. I VII We finally conclude in Sec. IVIII In this paper, 
the term "capsid" is used for both partial and complete 
spherical shells of viral proteins. The meaning should be 
clear from the context. 



II. THE ELASTIC MODEL OF HIV CAPSID 
BUDDING AND SELF-ASSEMBLY 

Let us consider a membrane-capsids system in which 
the concentration of Gags on the membrane, cq, is fixed. 
We assume all capsids assembled by Gags have the same 
average zenith angle a (see Fig. [2^) , and an average con- 
centration, n. n is related with a by the conservation of 
mass of Gags: 



n = c G - 



Mpta) 
A(a) 



CG- 



cos aa 



1 — cos a 



(3) 



where 



A{a) = 2irR 2 / smOdd = 2irR 2 {l - cos a) (4) 
Jo 

is the area of a capsid with zenith angle a, and Qg is the 
zenith angle of a single Gag (see Fig. [2Ji). Within this 
average description, it is convenient to think that the 
whole membrane surface is divided into identical cells, 
each contains a single capsid. The average size of these 
approximately circular cells, d, is given by the condition 



n{d/2) 2 n = 1. 



(5) 



Generically, the free energy density of the membrane- 
capsid system can be written as 



f = ne = n(e m + s c ), 



(6) 



where e is the free energy of one membrane cell. It in- 
cludes two parts: the elastic energy of the membrane, 
£ m , and the capsid energy e c coming from the Gag-Gag 
interaction and the Gag-membrane interaction. 

To calculate the elastic energy of the cell membrane, 
we use the standard Helfrich model 0, [TtJ where e m is 
the sum of two contributions from the bending energy 
and the stretching energy: 



dS 



' : (2H-C ) 2 + KG K 



dSa. 



(7) 



Here the integration with the area element dS is taken 
over the membrane surface, k and kq are the bend- 
ing rigidity and Gaussian bending rigidity, H and K 
are the mean and Gaussian curvatures, and Cq is the 



4 



spontaneous curvature of the membrane surface. Using 
the Gauss-Bonnet theorem, one can show that the total 
Gaussian curvature of the membrane surface is propor- 
tional to the total area of capsids, in the generic case 
when kq takes different values for the membrane at- 
tached to the capsid and the Gag-free membrane. Since 
the Gag concentration cq in our system is fixed, this term 
gives a constant in / and can be dropped from further 
consideration [l8j]. For a given Gag concentration cq, 
under our spherical capsid assumption, the shape and 
the total area of all capsids are fixed. Therefore the to- 
tal elastic energy of the membrane attached to capsids 
is also constant, and can also be dropped from consider- 
ation. As a result, the a-dependent contribution to e m 
comes from the integration over the Gag-free membrane 
surface only. In this region, we take the spontaneous 
curvature to be Co = 0, corresponding to normal lipid 
bilayer membranes. 

In consideration of the single capsid energy e c , since 
cq is constant, both the total Gag- membrane interaction 
energy and the bulk part of the Gag-Gag interaction en- 
ergy are constant. The only a-dependent contribution to 
e c comes from the rim energy of the capsid, due to the 
fact that the coordination number of Gags on the rim is 
not as many as Gags inside the capsid. Since the perime- 
ter of the capsid rim with zenith angle a is 27ri?sina, we 
set 



e c = r27r_Rsina. 



(8) 



The proportionality coefficient r can be considered as 
the "line tension" of the capsid. It is directly propor- 
tional to the strength of the Gag-Gag attraction and can 
be changed experimentally by mutations of Gags or by 
changing pH or salinity of the cell cytoplasm near the 
membrane. 

To proceed further, we take the "ideal capsids" ap- 
proximation when the distance between capsids is large 
and the membrane mediated interaction between them is 
negligible. Such an effective long-range interaction is pos- 



sible because the presence of the first capsid may change 
the deformation of the membrane around the second cap- 
sid and provides an effective interacting energy between 
the two. Qualitatively, this interaction is negligible when 
the capsid concentration n is small (the quantitative con- 
dition will be given in the next section) . Under this non- 
interacting capsids approximation, e m comes from the 
membrane deformation induced by a single capsid. 

The calculation procedure to find the free energy pro- 
file f(a) is as follows. We first minimize the membrane 
clastic energy e m with respect of all possible membrane 
shapes for any given capsid size a. Here it is convenient 
to use a cylindrical coordinate system (r, h, 4>) as shown 
in Fig. [2^ (the azimuthal angle <\> is not shown). With 
our assumption of (partial) spherical capsids, the mem- 
brane profile is independent on 0. As a result, one can 
use either the function h(r) or r(h) to parameterize the 
membrane. Correspondingly, the mean curvature and 
the area element can be written as Il9f 



H(r) 
dS 

or H(h) 
dS 



h'jrf + h'{r)+rh"{r) 

2rfl + h'(r) 2 ] 3 / 2 
ry/1 + h'(r) 2 drd(j); 
1 + r'{hf - r(h)r"(h) 
2r(h)[l + r'(h) 2 } 3 / 2 ' 

r{h)^\ + r'(h) 2 dhd<t), 



(9) 
(10) 

(11) 
(12) 



where h'(r) = dh/dr and h"(r) — d 2 h/dr 2 are the first 
and second derivatives of h with respect to r. Similarly, 
r'(h) = dr/dh and r"(h) = d 2 r/dh 2 are the first and sec- 
ond derivative of r with respect to h. Functionally min- 
imizing the membrane energy e m with respect to mem- 
brane shape r(h) or h(r), one obtains an elastic equation 
of the membrane shape, similar to the Euler-Lagrange 
equation derived from the least action principle in the 
classical mechanics. For the shape parametrization using 
r(h), Ss m /Sr = leads to the equation: 



r 12 o 14 o '6 /8 i // o '4 // o 16 II , r>„2 111 i 2 12 112 i o 2 14 1/2 r 3 f/3 

:[—r — Sr — 3r — r + rr — irr r — zrr r + 2r r — llr r r — 13r r r — or r 



2 r 2[ 1 + r /2]9/2i 

+30rVV' 3 + 4rVr( 3 ) + 8rW 3 ) + 4rW 3 ) - 20rVrM 3 ) - 20rVW 3 ) + 2rV 4 > + 4rVV 4 ) + 2rVV 4 >] 
1 + r' 2 - rr" 



[1 +r '213/2 



= 0, (13) 



where r^ = d 3 r/dh 3 and r^ = d 4 r/dh 4 are the third or 

and forth derivatives of r with respect of h. This equation , _ '(h)\ - cot (15) 

has to be solved together with the boundary conditions. ^ J I -R cos a — , )\r cos a ■ { ) 

On the rim of the partial spherical capsid, the membrane Far away from the capsid, the membrane becomes flat, 

itself and its slope must be continuous. We have we have 

h{r)\ Rsina = Rcosa, h'(r)\ Rsiaa = - tana, ; (14) h'(r)\ oo =0 (16) 
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r'(h)U 



oo. 



(17) 



Solving the elastic equation (|13p with the boundary 
conditions, Eq. (T5J) and {TTJl (or Eq. (O and Eq. ([I6j if 
/i(r) is used), one obtains the membrane shape that min- 
imizes e m . Substituting this shape into Eq. 0, one ob- 
tains the minimal e m (a). Putting its value into Eq. ©, 
one gets the total free energy density profile f(a). In 
general, the elastic equation, Eq. (fT3")) . is highly non- 
linear and numerical calculations are needed to obtain 
the exact membrane profile, as shown in Sec. [V] How- 
ever, in certain asymptotic limits, analytical solutions 
can be obtained which determine the scaling behavior of 
the system. This is done in the next two sections. 



III. ASYMPTOTIC SOLUTIONS OF THE 
MEMBRANE ELASTIC ENERGY 

In calculating the free energy profile, the most nontriv- 
ial part is to find the minimal e m (a), due to the nonlinear 
elastic equation involved. After the solution is found, it is 
straightforward to add the other part of the energy e c (a) 
and get /(a). Therefore we focus on the solution of min- 
imal e m in this section. Although not solvable in general, 
the problem do have analytical solutions in asymptotic 
limits. To a large extent, they determine the analytical 
behavior of the system, especially the scaling behavior of 
e m with the dimcnsionless parameter 



R 



(18) 



which characterize the relative strength of the surface 
tension to the bending rigidity. 



A. The small deformation solution 

A typical approach to consider the elastic deformation 
of the membrane is to take the small deformation ap- 
proximation which assumes |V/i| <C 1 [20| . Here we use 
the notation 



rd r + (p-drj, 



(19) 



in order to show similarity of the elastic equation to the 
linearized Poisson-Boltzmann equation later. Expanding 
with Vh and keeping terms of 0(V/i) 2 in 5e m = 0, we 
reach a linearized elastic equation which can be written 
as 



H 



V 2 H 



. 1, 
2 
H 



-Vh, 



0. 



(20) 



where we have introduced an important length scale in 
the problem, 



(21) 



It is the length scale beyond which the stretching en- 
ergy becomes more important than the bending energy. 
Notice that Eq. (|2"D|) takes exactly the same form as a 
linearized Poisson-Boltzmann equation in electrolytes or 
plasma 2l| . Therefore r s can be interpreted as an elastic 
screening length, similar to the Debye-Huckel screening 
radius. The local curvature H(r) induced by the capsid 
decreases when r increases and becomes exponentially 
small at projected distance larger than r s . This is a typ- 
ical linear solution of small deformation. 

Using boundary conditions Eqs. (TH| and (Til)]) , the spe- 
cial solution to Eq. (j2"0|) is given by 



h(r) 



R cos a + r s tan a 



K (r/r s ) - K (R sin a/r s ) 



Ki(Rsmu/r s ) 



h'(r) — —tana 
H(r) 



Ki{r/r s ) 



K\ (Rsma/r s ) ' 
tana Ko(r/r s ) 



2r s Ki (R sin a/r s ) 



(22) 
(23) 

(24) 



where K$ and K\ are the zero and first order modified 
Bessel function of the second kind. At r > r s , both 
K (r/r s ) and Kx(r/r s ) decay like y/r s /rexp(-r/r s ), 
and the deformation becomes exponentially small, as the 
meaning of r s suggested. 

Substituting this solution back to Eq. |7|), we get the 
minimal clastic energy of the membrane 



2 R sin a Kq [R sin a/r s ) 
E m = 7tk tan a — 



Ki (R sin a/r s ) 



(25) 



Notice that this energy is proportional to the dimcnsion- 
less parameter a = R^J a j k = R/r s . Here the inverse 
proportion to r s is again a generic feature shared with 
the theory of Debye-Hiickel linear screening [2l[ . 

The self-consistency of the small deformation approx- 
imation is warranted by |/i'(f)| < 1, or, according to 
Eq. (|23]) . | tan a | < 1. Therefore this solution is appli- 
cable in the whole range of r for a < tt/4 capsids only. 
On the other hand, at large distances far away enough 
from the capsid, the deformation of the membrane always 
becomes small enough such that the small deformation 
solution is applicable. In this sense, this solution can al- 
ways serve as a "far-capsid" solution for the membrane 
shape, although the formula for e m in Eq. (|25|) is not valid 
in general. It describes the universal decaying behavior 
of the deformation when the deformation itself becomes 
small enough. We can formally define a characteristic 
distance r c through 



\h'(r c )\ = l, 



(26) 
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beyond which the small deformation solution is valid. r c 
will be useful later when we discuss the complete solution 
to the problem. 

With the small deformation solution in hand, we are 
now ready to derive a quantitative condition for the 
ideal capsid approximation introduced in the last sec- 
tion. Clearly, when the average projected distance be- 
tween capsids, do, is much larger than r s + 2R, the mem- 
brane mediated interaction between capsids is negligible, 
since the deformations of the membrane by the capsids 
at distance larger than r s are screened out. In this case, 
most of the membrane surface is flat, so do — d (notice 
that d is measured along the membrane surface which in 
general is larger than do measured along r axis). Thus 
according to Eqs. (J3j> and ([5|), the ideal capsids approxi- 
mation is valid when 



do 



2sin(a/2) 



r s + 2R (r s + 2i?)^r7j^sin(a G /2) 



> 1. 



(27) 



In this work, we assume cq is small enough and this is 
always the case. 



B. The catenoid solution 
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FIG. 4: (Color online). In the limit of small surface tension, 
the optimal membrane shapes (thin line, blue online) around 
the capsid (thick line, red online) are catenoids, as shown for 
different capsid sizes. 



When the surface tension a = 0, or, r s — > oo, again an 
analytical solution is available 22J. In this case, the sec- 
ond integral in Eq. (|7|) is zero. Our problem of finding the 
minimal e m is reduced to a minimal surface problem in 
differential geometry [l9(. Namely, we look for the solu- 
tion to the equation H = [23( . The only solution under 
the rotational symmetry of our problem is the catenoid 
solution, first discovered by Eulcr in 1740 (24|. 

In this case, due to the possible multiple values of h 
at the same r, it is better to use the r(h) representation. 
H is then given by Eq. (TPTj) . Using boundary conditions 
(fl"5|) and (fl"7|) . the special solution to H = is 



(28) 



"(h) = Rsin 2 a 



cosh 



h — R cos a — R sin 2 a arcsinh(cot a) 
R sin 2 a 



The catenoid shapes for various a are depicted in Fig. [4] 
In this catenoid shape, the elastic energy e m achieves its 
absolute minimum, zero. 

The catenoid solution is a solution to a nonlinear differ- 
ential equation. It involves large deformations which can 
not be characterized by the linear solution discussed in 
the last subsection. Although exact only when r B — > oo, 
this solution is still useful for large but finite r s [2f|. In 
fact, since there are no other length scales in the elastic 
equation (| 13(1 (i? only shows up in the boundary con- 
ditions), a large r s actually means r s ^> r. Therefore 
in the region of r -C r s , the catenoid solution should 
work asymptotically. In this sense, this solution can al- 
ways serve as a "near-capsid" solution for the membrane 
shape, although e m = is not true in general. The char- 
acteristic length beyond which it fails is simply r s . 



In case of r s — > oo and a <C 1, both the catenoid 
solution and the small deformation solution work in all 
range of r. Indeed they become identical. 



C. Membrane elastic energy at two asymptotic 
limits 

The two solutions discussed in the last two subsections 
determine the analytical behavior of the system to a large 
extent. When ? « 1, they can be combined to get the 
analytical expression of the minimal e m . In general, they 
determine the scaling behavior of e m with respect of a. 
We separate our discussion into two opposite limits of 
small and large a, which can be called the soft membrane 
regime and the stiff membrane regime. Here "soft" means 
easy to stretch, "stiff" means the opposite. 

In the soft membrane regime, a <^ 1 or R <C r s , the 
catenoid solution is valid near the capsid when r -C r s . 
Calculating r c using Eqs. ([28]) and (|26|) . we get 



V2R 



■ 2 

sm a. 



(29) 



We see r c <C r s . Therefore the valid regions of the two 
asymptotic solutions (one is r > r c , the other is r < r s ) 
overlap largely and we can combine them to get a com- 
plete solution to the optimal membrane shape. Quanti- 
tatively, we artificially choose a projected distance some- 
where between r c and r s , say ^/r c r s . For r < ^r c r s , 
the catenoid solution is used. For r > ^/r c r s , the small 
deformation solution is used. Notice that the special so- 
lution of the small deformation has to be calculated using 
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the continuity conditions for h(r) and h'{r) at y/r c r Sl de- 
rived from the catenoid solution. As a result, we have an 
analytical expression for the optimal membrane shape 
continuously from the edge of the capsid to infinity. The 
corresponding e m , keeping the leading order terms in the 
small parameter a, is given by 



■ 4 R\ r s 
7TKsm a— tt In — . 

r 2 R 



(30) 



A. The soft membrane regime 

In the soft membrane regime, a = R/r s -C 1. Substi- 
tuting Eqs. © and (|25J) to Eq. ©, we have 



/ = kc g (1 - cosa G )f 



kcq(1 — cosctc)TT cot t|-(2t + ct 2 In i sin 3 a), 



(34) 



When a <C 1, this result agrees with the small deforma- 
tion solution in Eq. (|25[) in the same regime of small a. 
In this limit, we have e oc <t 2 ln(l/(r). 

In the stiff membrane regime, a 1 or iJ > r s . 
Since r ^> r s always, the "near capsid" region where the 
catenoid solution holds disappears. On the other hand, 
for a < 7r/4, the small deformation solution is valid in 
the whole range of r. The membrane elastic energy is 
given by (|25p . which in this limit reads, 



where we have introduced the dimensionless free energy 
density / for convenience, f(a) is plotted schematically 
in Fig. [3] When r is large, the only minimum of the free 
energy density is at a — > tt (the left panel of Fig. [3]). On 
the other hand, when r < 0.2er 2 ln(l/cr), a local minimum 
at the capsid size, ao, appears ( the right panel of Fig. [3]). 
Correspondingly, the threshold line tension at which the 
local minimum in the free energy density appears is 



0.2i?erln- 



(35) 



2 • R 

£ m =7T«;tan asina — . 



For a > 7r/4 capsids, a rough estimate of r c using Eq. 
gives 



R sin a 



In I tan a I 



(31) 

Since transcendental equations are involved in mini- 
mization of /, it is not easy to get the analytical expres- 
sion about this local minimum in general. However, ag 
and the corresponding /o can be estimated in a linear 
approximation. Assuming ao is achieved at small a, we 
can expand / and keep only the leading order terms in 
(32) a . We get 



Since R 3> r s , for most of a, we expect that the small 
deformation solution starts to work at places close to the 
capsid. Probably because of this, the scaling behavior of 
e m ex a is preserved even at large a, as shown by the 
numerical result (see Sec. PVT) . 



/ = h 27R7 2 ln-a 3 . 

a a 



Taking df /da = 0, we have 



a 



o- 2 ln(l/a) V/ i?CTln(V^/Vi?) 



(36) 



(37) 



IV. ANALYTICAL RESULT OF THE TOTAL 
FREE ENERGY DENSITY 



After the information about the minimal s rn (a) is 
known, we can add the line tension energy e c (a) to it and 
consider the total free energy density f(a). The presence 
of e c introduces the second dimensionless parameter to 
the problem, 



The fact that this is a minimum rather than a maximum 
is confirmed by d 2 f /da 2 \ ao > 0. For r < t c at which ao 
shows up, this result is indeed much smaller than one, 
consistent with the initial assumption that «o < 1. In 
the same limit, 



fo 4-7T y/r 2 a 2 ln(l/a) = 47r 




a In ■ 



R 



(38) 



which characterize the relative strength of the line ten- 
sion on the capsid rim. In this section, we derive sev- 
eral simple scaling behaviors of the system, depending 
on the two dimensionless parameters a and r. We again 
separate our discussion into the soft and stiff membrane 
regimes corresponding to small and large a. 



B. The stiff membrane regime 

In this case, we do not know the form of the membrane 
elastic energy e m for large a. Still, in the same spirit of 
linear analysis, we can assume that there is a minimum 
of / at small a, and use the small deformation solution 
Eq. (|3ip for s m . Notice that the minimum found in this 
way is only a local minimum, since we did not include 
the information of large a. 
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As a result, we have 
f 



Li . o \ 

7r cot — (2r + a tan a) 



Attt 
a 

Taking df /da = 0, we get 



a 



2naa. 



(39) 



(40) 

v fto 

It is a minimum since d 2 f /da 2 \ aa > 0. For this result to 
be meaningful, r <C cr must hold, which will be checked in 
comparison with the numerical result. The correspond- 
ing free energy density is 




fo = 4-7TV 2(JT = AuR 



A(TT 2 



(41) 



V. NUMERICAL RESULT AND DISCUSSION 

In order to verify our analytical understanding and get 
the complete solution to the problem, we solve the non- 
linear elastic equation derived from Se m = numerically. 
Our computation procedure follows Refs. jf| 1261 ]. This 
numerical solution is then combined with s c to give the 
total free energy density /. In this section, we show the 
numerical result, compare it with the analytical formulas, 
and discuss the meaning of our results. 



In 



-2 



2 1 < 



v 



I -Li o 



U -i o 1 » 




0.5 1.0 1.5 2.0 2.5 3.0 



a 



is that the elastic energy profile always takes a "sand 
dune" shape, where two minimums, zeros, are achieved 
at a — > 0,7r, and a maximum shows up in the middle of 
a. Physically, this energy profile comes from the need of 
matching boundary conditions at the edge of the capsid 
and at infinity. The membrane deformed by the capsid 
edge at one end has to become flat far away from the 
capsid. At a — > and a — > 7T, the membrane is not 
deformed at all, and the elastic energy is zero [13] ■ While 
for a close ton/2, the membrane is almost vertical at the 
edge of the capsid, and a large amount of elastic energy 
is needed to bend it flat. 

Secondly, we see clearly two kinds of asymptotic be- 
haviors of e m depending on the parameter a = R/r s . In 
the stiff membrane regime, cf ^> 1, the energy is propor- 
tional to a as shown by the collapse of the data points 
to a single curve with a varying from 10 2 to 10 5 . The 
maximum of the energy is achieved at a m ~ 0.7n. a m is 
a nonlinear result and can not be calculated analytically. 
However, the proportionality of e m to a is a small defor- 
mation result as shown in Eq. (|31|) . In the soft membrane 
regime, cr <C 1, the energy is proportional to a 2 hi(l/a), 
shown again by the collapse of the data points with a 
varying from 10~ 2 to 10~ 5 . Here the collapse is not as 
pronouncing as in the other regime mostly due to the 
larger numerical error in dealing with smaller e m . The 
absolute value of e m in this regime is smaller at least in 
four order of magnitude than in the other regime. The 
maximum of the energy here is arrived at a m — n/2 
and the curve becomes symmetric about a m . These fea- 
tures agree with our small a solution originating from 
the catenoid solution. In fact, Eq. (|30|) fits the numerical 
data reasonably well, with an additional factor 1.3. 



£ in 
cf 




FIG. 5: (Color online) Numerical result of the dimensionless 
membrane elastic energy e m — £ m /s as a function of a. The 
eleven sets of data points are at a = 10~ 5 , 10" 4 , 10" 3 , 10 5 . 
They are labelled correspondingly as —5, —4, —3, ...5. The left 
axis em/a is for all a > 1 data points (red on line), while the 
right axis Em/cr 2 ln(l/a) is for all a < 1 data points (blue on 
line), as indicated by the two arrows. The curve represents the 
analytical asymptotic solution (|30p with an additional factor 
1.3 (blue on line), fitting the data points for 5 < 1. 

The direct numerical result of e m is plotted in Fig. 03 
where for convenience we used the dimensionless elastic 
energy e m = e m /n. The first important thing to notice 



FIG. 6: An effective "phase diagram" in the plane of two di- 
mensionless parameters, a — R/r 3 and f = Rt/k. In the 
upper-left part, the free energy density decreases monoton- 
ically with a, while in the lower-right part, it has a local 
minimum, as shown in the insets. The numerical data points 
mark the "phase boundary" at which the local minimum ap- 
pears. The dotdashed line and the solid lines fit the data 
points using r = O.IIct and r = 0.065<r 2 ln(l/5) respectively. 

The scaling of e m with cr suggests a simple way to 
do the numerical calculation to the free energy density. 



9 



When cr > 1, 

7=e c 



e m = 2ir cot — r + 31 (a)er 



(X T 

27rcot-- +31 (a) 
A a 



(42) 



where e c = £ c /k and 31(a) is some function given by the 
numerical computation. According to the last equality, 
up to an overall constant a, f is completely determined 
by only one parameter t/ct. Similarly, when a -C 1, 



5 2 ln4 
a 



a 

2n cot — 



2 a 2 ln(l/a) 



32(a) 



(43) 



where 32 (a) is again given by numerical computation, al- 
though we know it from our analytical result Eq. (|34|) . 
In this regime, / is determined by one parameter 
t/ct 2 ln(l/5). Below in studying the local minimum of 
/, we therefore consider a single parameter dependence. 



a 



<x»l 




a 



<T«1 



0.00 0.01 0.02 0.03 0.04 0.05 

f/<T 2 ln(l/<r) 



FIG. 7: The capsid size ao as a free energy local minimum is 
shown as a function of r/a and t/ct 2 ln(l/er) at two limits of a. 
The dots are numerical results taken at a = 10 2 , 10 3 , 10 4 , 10 5 
for the upper panel and a = 10~ 2 , 10~ 3 , 10~ 4 , 10~ 5 for the 
lower panel. The curves are analytical results of Eq. (|40|) 
at 5 > 1 and Eq. (|37[l at 5 < 1 with additional numerical 
factors of 2 and 1.5 respectively. The range of f plotted cor- 
responds to the lower-right "phase" in the "phase diagram" 
of Fig. © 

For all a and r, we get two different types of free en- 
ergy density profiles as shown in Fig. [31 consistent with 



the analytical result for small a. The global minimum of 
the free energy density is always at a — > 7r. Physically, 
the line tension energy prefers the shortest length of the 
capsid rim, which is zero for complete capsids (cv — ► tt). 
When t is very large, the line tension energy dominates, 
and the free energy density / decreases with a monoton- 
ically to zero, as shown in the left panel of Fig. [3j On 
the other hand, when r is small, due to the maximum 
of the membrane elastic energy e m , a local minimum at 
the capsid size, ao, shows up in the free energy den- 
sity, as shown in the right panel of Fig. [3j It is useful 
to draw a "phase diagram" on the plane of a and r as 
Fig.[S]to show this qualitative difference in the free energy 
density profile. The lower right region of Fig. [6] corre- 
sponds to value of the parameters (a, r) where capsid 
budding can be kinetically trapped. The two lines fit the 
"phase boundary" at large and small a with r = 0.1 la 
and t — 0.065<7 2 ln(l/cr) respectively. As one can see, 
there is a very good agreement between numberical re- 
sults and our scaling formulas for a in two asymptotic 
limits. According to the numerical fits, the threshold r 
at which the local minimum in the free energy density 
shows up are 



T C = 0.11V 

when 5> 1, and 

r r = 0.065ifo-hv 



(44) 



(45) 



when a <C T The later formula agrees with our analytical 
result Eq. ()35[) with a numerical factor 3 difference. 

The possible local minimum in / (the right panel in 
Fig. [3j suggests that budding may be trapped kinetically 
at the capsid size ao- Numerical and analytical results of 
ao are shown in Fig [71 The analytical curves are drawn 
using Eq. (gDJ> at a > 1 and Eq. (J27J at a < 1, with addi- 
tional numerical factors of 2 and 1.5 respectively. There 
is a large deviation between analytical and numerical re- 
sults at large r. This is the parameter regime where the 
linear approximation is no longer valid. 

The kinetic trapping becomes significant if the barrier 
in the free energy density is large. In Fig. [51 numerical 
and analytical results about this barrier are plotted. For 
the local minimum /o, up to a order one numerical fac- 
tor (1.7 and 3.2), our analytical expressions Eq. (|4"Tj) at 
a ^> 1 and Eq. (|38|) at a <C 1 remains reasonable approx- 
imation. We can not estimate the maximum f m , which 
is in the nonlinear regime. However, in the most impor- 
tant regime of small r and large barrier, Fig. [8J shows 
that the main contribution to f m comes from the mem- 
brane elastic energy e m (the value of f m at t = 0). In 
this regime, the additional contribution to f m from the 
line tension energy e c is negligible and f m is almost a 
constant. Combining the numerical result of f m and the 
analytical result of /o with proper numerical factors, we 
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f/(T 2 ln(l/(T) 



FIG. 8: (Color online). The local minimum fo, maximum f m , 
and barrier f m — fo of the free energy density. The values of 
a and f plotted are the same as in Fig. [7] The circles (blue 
online) are the numerical result of fo, fitted by the solid lines 
(blue online) using Eq. (|4"Tj) at a > 1 and Eq. (O at a <C 
1 with additional numerical factor 1.7 and 3.2 respectively. 
The squares (green online) are the numerical result of f m , 
marked by the dashed lines (green online) at their zero t 
values. The triangles (red online) are numerical result of the 
barriers, f m — fo, fitted by the dotdashed lines (red online) 
using Eqs. (gg} and (gTJ . 



get the asymptotic formulas for the barrier at t < 1, 



fm — fo 



7.8a - 6.87rV2CTr 



7.sIlJ--6.8irRj— r , (a > 1). (46) 



f m -f ~ 4.4a 2 ln(l/a) - 12.8tt ^a 2 ln(l/oF) 



= 4.4i? 2 -ln^ 12.87r-Wr 2 crln- V 



K 

(a « 1). 



R 



R 



(47) 



VI. KINETICS OF BUDDING AND PARTIAL 
BUDDING 

As discussed in the previous sections, with a finite Gag- 
Gag attraction, budding always proceeds to completion 
thermodynamically. However, when this attraction is 
weak, or r is small, a metastable state of partial bud- 
ding appears at a smaller capsid size ao (see Fig. [3]). It is 
therefore possible that the budding process is kinetically 
trapped at ao- In this section, we discuss this kinetic ef- 
fect and make connections of our theory to experiments. 

Let us first estimate the values of parameters. A 
normal plasma membrane has k ~ 20 — 40keT and 
a — 0.5 — 2pN/nm = 0.12 - 0.48k B T/nm 2 %2§L On the 
other hand, typical HIV have R ~ 60 — 80nm PJ. Conse- 
quently, a — Ry/a/n ~ 10 and only the stiff membrane 
regime with a ^> 1 is relevant for HIV. Actually in the 
opposite regime of a <C 1, one can show that the typi- 
cal energy scale e m ~ na 2 ln(l/er) is comparable to keT 
and thus not important in the room temperature. In this 
section, we therefore focus on the stiff membrane regime 
only. 

In order to see if budding can be kinetically trapped at 
the local minimum ao (see Fig. [3]), we have to study the 
budding kinetics and calculate the kinetic barrier. For 
this purpose, let us employ the standard kinetic picture 
of the first order phase transition [29|, [3(|, correspond- 
ing to the transition from a free-Gags phase to an aggre- 
gated Gags phase where Gags self-assemble into complete 
viral capsids. At the initial stage of aggregation, the con- 
centration of free Gags is large, Gags coagulate to form 
dimers. Dimers coagulate with free Gags or other dimers 
to form larger Gag clusters (small capsids). This initial 
coagulation or nucleation is a fast process and is not a 
rate limiting step in retroviral budding. Soon free Gags 
are significantly depleted, and the main kinetic pathway 
for growth of capsids is for them to diffuse and merge 
with each other. We will be concern with this later stage 
of coagulation. For simplicity, we work with the domi- 
nant capsid size, a(t) [with concentration n(t)], assuming 
these typical capsids carry all the mass of membrane- 
bound Gag proteins. 

Let us start with the case when a(t) is still small so 
that the energy barrier for merging of capsids is smaller 
than kpT. This is the regime of the well known diffu- 
sion limited aggregation ,31|. The rate of the capsid area 
A(a, t) incretion is proportional to the probability that 
two capsids diffuse and merge with each other. The ki- 
netic rate equation reads 



dA(a, t) 
Jt 



2?ri? sin a(t)A{t)DVn{t)\ R sina , (48) 



The largest barriers are achieved at r = or fo = 0. 



where D ~ kgT '/rjRsm.a(t) is the lateral diffusion con- 
stant of a capsid on the membrane and Vn(t)|# S i na is the 
gradient of the concentration n(t) on the edge of the cap- 
sid. This gradient can be estimated assuming a steady 
state in the diffusion and taking the adsorbing boundary 
condition at the edge of the capsid and a given capsid 
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concentration [Eq. ([3])] far away from the capsid. Solving 
the diffusion equation with these boundary conditions, 
we find 



Vn(t)|ij si na 



c G A(a G ) 
A{a) i?sina 



(49) 



Substituting these relations and Eq. (j4|) into Eq. (|48p , we 

get 



a(t) sin2a(<) t 



T drff 



2 



sin 2a G 



where 



Tdiff = f]R 3 /Tc G A(a G ) 



(50) 



(51) 



is the time scale of diffusion proportional to the viscosity 
?y of the membrane. In the small a regime corresponding 
to a small kinetic barrier, this equation can be written as 



a(t) = {Zt/ Tdiff 



~3 \l/3 
a G) i 



(52) 



which is a slow function of time. 

The regime of diffusion limited growth stops when 
the kinetic barrier between approaching partial capsid is 
much larger than k B T . At a later time, a different growth 
regime of Lifshitz-Slezov (LS) comes into play [23]. In 
this mechanism, the growth is no longer due to colli- 
sion and merging of partially budded capsids. Instead, 
smaller capsids shrink and release individual Gags. These 
Gag are absorbed into larger capsids, leading to their 
growth. This process of releasing and adsorbing of indi- 
vidual Gags (the so-called coalescence) has much smaller 
kinetic barrier than the barrier to capsid merging in this 
later stage. The growth of capsid size in LS regime is the 
same as that of diffusion limited growth [29j. However, 
the rate constant tls depends exponentially on the acti- 
vation energy to release individual Gag proteins from a 
capsid 



t ls oc exp(-\e\/k B T) 



(53) 



where e is the binding energy of Gag in a capsid, which 
itself is also a function of the Gag-Gag interaction. 

The kinetic picture described above is good when r > 
t c and the free energy density decreases monotonically 
with increasing a (the left panel of Fig. [3]). However, 
when t < t c and a local minimum ao appears in the 
free energy density (the right panel of Fig. [3]), the above 
picture must be modified. For the cluster growth, either 
in the diffusion-limited regime or in the LS regime, the 
growth of the cluster size always reduces the free energy 
of the system. On the other hand, for the capsid growth 
of retroviral budding, after the capsid size reaches ao, 
the system free energy increases when the capsids grow 
further. For a > ao, the growth of capsids is determined 
by the ability to tunnel through the kinetic barrier related 
to f m — fo ( see Fig. UJ. The detail analysis of the rate of 
capsid growth for a > ao is a very interesting problem 



by itself, requiring understanding of membrane energetics 
when a partially budded capsid absorbs other capsids or 
many individual Gags to increase its size from ao to a m . 
These calculations are beyond the scope of this paper 
and we will leave the detail treament of capsid growth in 
this case to a future study. Nevertheless, one can expect 
the rate of such process to inversely proportional to the 
exponential of the energy barrier 



e i oc exp[-(/ m - f )/nk B T], 



(54) 



where (/ m — fa)/n is the energy barrier of a membrane 
cell with a single capsid in it. According to Eqs. ([46]) 
(see also the upper panel of Fig. [8]), the maximum energy 
barrier is achieved at fo = or f = 0. Using Eqs. ([3J 
and (1341) . it can be written as 



Em. — 



In 



= k(1 - cosa m )/„ 



(55) 



where f m is given by Eqs. (|46p . and a m is the correspond- 
ing capsid size. A more useful expression of E m can be 
got if one recognizes that E m is nothing but the maxi- 
mum of E m shown in Fig. [5] Using the numerical result 
of that figure, we get 



En 



11.5kcj = 11.5i?y / Ko r . 



(56) 



Clearly, E m ^> keT for a > 1. For example, for 
R = 70nm, a = 0.24k B T/nm 2 and n = 20k B T, we get 
E m = 1765keT. The true energy barrier is smaller than 
E m since f > 0. In the regime of small r, according to 
Eq. |@BJ), it is 



E ~ k(1 - cos Om)(f m - fo) 




In experiments, if one treats k and a as constants then 
according to Eq. ([56")) the larger the retrovirus size R, 
the bigger the kinetic barrier. On the other hand, the 
line tension r is directly proportional to the strength of 
the Gag-Gag attraction and is experimentally adjustable 
through mutation of the late domain on the Gag pro- 
tein, binding of other molecules to Gags or changing the 
pH, salinity of water solution near the membrane [2lll3l|. 
As we know, the closest approach distance between two 
Gag proteins is about lOnm jl]. If Gags are densely 
packed on the capsid, this gives r ~ 0.5kBT/nm ~ 2pN 
for normal retroviral capsids. Theoretically, in order 
to have a local minimum in the free energy density 
and trap retrovirus budding kinetically, we must have 
t < t c (see Fig. [3]). For a normal cell membrane with 
k = 20k B T and a = 0.24k B T/nm 2 , using Eq. (JUJ), 
r c ~ 0.24kBT/nm = lpN. Therefore for normal capsids, 
r > t c , and budding easily proceeds to completion (see 
the left panel of Fig. [3]). One the other hand, r is bigger 
than t c only by a factor of 2. Therefore HIV budding 
can be fairly easily trapped at a partially budded state 
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with capsid size ao by reducing the Gag-Gag interac- 
tion strength such as mutation of a single domain on the 
Gag protein. The kinetic barrier E appeared at a — ao 
can be much larger than keT, and the time scale for 
capsid growth beyond ao, Ttunneh is exponentially large. 
Qualitatively, this trend is consistent with experiments 
on mutation of the late domain of Gag proteins @, 0]. 
Numerically, we know that ao < 0.9 ~ 0.37T (see the 
upper panel of Fig. (T)) . It agrees with experiments rea- 
sonably well, although there are many additional factors 
that we neglected in our treatment such as local varia- 
tion in membrane elasticity due to raft structures or the 
presence of other proteins in in-vivo assembly and bud- 
ding. More controlled experiments are needed to to verify 
the dependence on the membrane rigidities and Gag-Gag 
attraction of ao given by Eq. (HO)) . 

VII. CONCLUSION 

In this paper, we developed an clastic model of HIV 
(and retroviruses in general) budding and self-assembly 
on the elastic membrane. We studied the free energy 
profile of the system as a function of the capsid size 
a. We showed that although always thermodynamically 
favorable, complete budding and assembly may not be 
achieved if the Gag-Gag attraction is weak. In practice, 
for normal biological conditions, the Gag-Gag attraction 
is strong enough and HIV budding and assembly always 
proceed to completion, as it should be. On the other 
hand, it is fairly easy to trap HIV budding to a par- 
tially budded state with capsid size ao by reducing the 
Gag-Gag attraction. This can be done through the mu- 
tation of late domain on the Gag protein or binding of 
other molecules to Gag. In principle, the trapping is also 



possible by increasing the membrane rigidities, although 
this is not easy to do in vivo. Our theory agrees with rea- 
sonably well with experimental results. However, experi- 
ments with better controlled environments are needed to 
verify various aspect of the theory. 

The most interesting point of our model is probably 
that it provides a unique self-assembly mechanism. Not 
like self-assembly of other viruses or colloids, HIV assem- 
ble and bud concurrently on the membrane. Therefore 
the membrane elastic energy plays an important role in 
the assembly process. For example, the kinetic barrier 
which traps the HIV budding essentially comes from the 
membrane elastic energy around the capsids. In fact, 
our model developed for HIV budding and assembly can 
be very well applied to other situations. For example, 
for a given concentration of membrane-bounded proteins 
with a fixed spontaneous curvature, this kind of bud- 
ding and assembly phenomenon should also exist and 
can be explained using our model. In this situation, it 
may be easier to change the membrane properties and 
protein-protein attraction in vitro to verify our theory 
more quantitatively. Due to the interplay between the 
membrane elastic energy and the Gag-Gag attraction en- 
ergy, the kinetics of retrovirus budding is an interesting 
problem by itself, as discussed in Sec. EB We plan to 
address this question in more detail the near future. 
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